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ABSTRACT 

The instability of a compressible flow past a wedge is investigated in the hypersonic 
limit. Particular attention is given to Tollmien-Schlichting waves governed by triple-deck 
theory though some discussion of inviscid modes is given. It is shown that the attached 
shock has a significant effect on the growth rates of Tollmien-Schlichting waves. Moreover, 
the presence of the shock allows for more than one unstable Tollmien-Schlichting wave. 
Indeed an infinite discrete spectrum of unstable waves is induced by the shock, but these 
modes are unstable over relatively small but high frequency ranges. The shock is shown to 
have little effect on the inviscid modes considered by previous authors and an asymptotic 
description of inviscid modes in the hypersonic limit is given. 
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1. INTRODUCTION 


Our concern is with the instability of hypersonic flow around a wedge of small an- 
gle. The work is motivated by recent interest in the development of hypersonic vehicles. 
We do not account for real gas effects even though some of the vehicles in question will 
certainly operate at speeds where such effects are important. Here we shall concentrate 
on Tollmien-Schlichting waves governed by triple-deck theory though some discussion of 
inviscid modes will be given. In a related calculation, the Gortler vortex instability mech- 
anism at hypersonic speeds is considered. There it is shown that centrifugal instabilities 
have their structure simplified considerably in the hypersonic limit. Here we show that 
this is not the case for Tollmien-Schlichting waves though inviscid modes have a relatively 
simple structure. In fact the logarithmically small layer at the edge of the boundary layer 
which controls inviscid modes is precisely the same layer where Gortler vortices become 
trapped at hypersonic speeds. It would therefore seem that the nonlinear interaction be- 
tween inviscid and centrifugal instabilities is an important problem to be considered in the 
hypersonic limit. 

Before giving more details of the problem to be considered here we shall give a brief 
review of some relevant previous calculations on compressible stability problems. For an 
excellent review of viscous and inviscid stability properties of compressible boundary lay- 
ers the reader is referred to the article by Mack (1987). Perhaps the main feature of 
compressible boundary layers of practical importance is that they are unstable to both 
inviscid and viscous instability waves. The inviscid modes have wavelengths comparable 
with the boundary layer thickness whilst the Tollmien-Schlichting waves have longer wave- 
lengths. Either the successive approximation procedure of G aster (1974) or the formal 
asymptotic description based on triple-deck theory, as used for incompressible flows by 
Smith (1979a, b), can be extended to the compressible viscous instability problem in a rou- 
tine manner. Such a calculation was given by Gaponov (1981) using essentially Gaster’s 
approach whilst more recently Smith (1988) has applied triple-deck theory to the viscous 
modes of compressible boundary layers. A significant result of Smith (1988) is that when 
the free stream Mach number is increased there is critical size of Mach number in terms of 
the (large) Reynolds number at which the Tollmien-Schlichting downstream development 
takes place on a lengthscale comparable with that over which the basic state develops. At 
this stage the waves cannot be described by any quasi-parallel theory and the evolution 
can only be described by numerical integration of the governing linear partial differential 
equations. This is precisely the situation in the incompressible Gortler stability problem, 
see Hall (1983); perversely in the hypersonic limit the most significant Gortler vortices lose 
this property and become only weakly dependent on nonparallel effects. A question of some 
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importance raised by Smith’s work is whether the failure of the quasi-parallel approach at 
high Mach numbers means that many of the parallel flow calculations in this regime at 
finite Reynolds numbers are in error. Because Smith’s prediction is based on a double high 
Reynolds number and Mach number limit, the regime at which this failure occurs must be 
identified by numerical integration of the governing linear partial differential equations. 

The inviscid modes of instability of a compressible boundary layer have been well- 
documented by Mack (1987). In fact, there can be unstable two- and three-dimensional 
modes and neutral modes associated with a generalized inflection point and non-inflectional 
neutral modes. At Mach numbers above three it is the inviscid modes which have the high- 
est growth rates and therefore presumably dominate the transition problem in compressible 
boundary layers. However, previous calculations for viscous and inviscid modes have taken 
little account of the presence of shocks in the flow field (but see Petrov 1984) , though at 
high Mach numbers there is no question that they play an important role in determining 
the basic state. A primary aim of the present calculation is to gain some insight into 
the role of shock waves in the formation of viscous and inviscid modes in a compressible 
boundary layer. 

The particular configuration which we investigate is the hypersonic flow around an 
aligned wedge of semi-angle 0. The inviscid flow in this case is a uniform state either side 
of straight shocks which make an angle <j> = a — 0 with the wedge. At the wedge a boundary 
layer is set up and the fluid velocity is reduced to zero inside this layer. We restrict our 
attention to the case when the wedge is insulated though our calculations can readily be 
extended to the isothermal case. 

In order that Tollmien-Schlichting wave disturbances to this flow can be treated in a 
quasi-parallel manner we make the Newtonian approximation and take the distance of the 
shock from the wall to be comparable with the upper deck scale over which Tollmien- 
Schlichting waves governed by triple-deck theory will decay. At the shock we derive 
linearized boundary conditions which the Tollmien-Schlichting wave must satisfy. This 
condition effectively changes the eigenrelation from that discussed in the absence of shocks 
by Smith (1988). We show that the shock has a significant effect on the growth rate of the 
mode corresponding to that discussed by Smith. In addition we show that the presence 
of the shock allows for the existence of much more amplified modes trapped between the 
wedge and the shock. These modes occur at relatively high frequencies but are unstable 
over very short ranges of frequency. 

Some discussion of inviscid disturbances is also given. At high Mach numbers we give 
the appropriate asymptotic structure of the so-called acoustic modes, and show that they 
will be influenced little by the presence of the shock. The results which we give in the 
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hypersonic limit are in remarkably good agreement with Mack’s (1987) results even at 
relatively low Mach numbers. 

The procedure adopted in the rest of this paper is as follows: in §2 we formulate the 
appropriate sealings for Tollmien-Schlichting instabilities of the hypersonic flow past a 
wedge. In §3 the triple-deck equations for such modes are derived and in §4 the dispersion 
relationship appropriate to the linearized form of these equations is given and discussed. 
In §5 and draw some conclusions about the results of §4 and show how the appropriate 
inviscid stability modes develop in a hypersonic boundary layer on a wedge. Finally in 
an Appendix we derive the shock conditions appropriate to a disturbance with arbitrary 
length scales and use them to find the simplified shock condition for a Tollmien-Schlichting 
wave. 


2. FORMULATION 


The basic flow whose stability we examine is illustrated in Figure 1. For simplicity, 
the wedge is taken to be symmetrically aligned with an oncoming supersonic flow with 

A 

velocity magnitude U. Shocks of semi-angle a develop from the tip of the wedge (the 
acute angle between the shock and the wedge is <t> = o — 0). Quantities upstream of the 
shock are indicated by the subscript u , and quantities in the so-called ‘shock-layer’ between 
the shock and wedge by the subscript s. Cartesian axes are introduced with the x and y 
coordinates aligned and normal with the upper surface of the wedge, and the z coordinate 
in the spanwise direction. The corresponding velocity components are u = (u,v,u>), while 

A A A 

t, p, p, T and h are used to denote time, density, pressure, temperature and enthalpy 
respectively. We will assume that the fluid is a perfect gas with ratio of specific heats 7; 
then the upstream Mach number, Af u , is given by 


where the sound speed a u is given by 



a u 


(2.1a) 


*\ = ^ = ft - 1)K. ( 2.1b ) 

Pu 

The inviscid solution for this flow configuration consists of uniform quantities on either 
sides of straight shocks (e.g. see Hayes and Probstein, 1966). Specifically, 


e = 


A 

Pu 

A 

Ps 



2 

(7 - 1 )Ml sin 2 o ) ’ 


(2.2a) 


f^- = 1 + 7 Ml sin 2 a(l - e), 


(2.2 b) 
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l + ^(7-l)(l-e 2 )M u 2 sin 2 a, 


(2.2c) 


tan <f> = e tan a, (2.2 d) 

(«||)u= (tty), = Ucosa, (ui) u = — U sincr, («x)* = — e 17 sin a, (2.2e) 

where tty and u± are the velocity components parallel and perpendicular to the shock. From 
(2.1) it follows that between the shock and the wedge the fluid velocity has magnitude 


U, = U cos a(l + e 2 tan 2 o) 5, 


(2.3) 


and the local Mach number M, is given by 

2 _ M l cos2 <7(1 + £ 2 tan 2 a) 

* “ H-|(7 -l)(l-e 2 )M u 2sin 2 a’ 

or equivalently 


(2.4a) 


Ml = 


Ml 


cos 2 o(l + e 1 tan 2 a) — £(7 — 1)(1 — e 2 )M 2 sin 2 a 


(2.46) 


Note that (2.3) specifies the slip velocity along the wedge, which in viscous flow necessitates 
the existence of a boundary layer. 

Before examining this boundary layer in detail we introduce a non-dimensionalization 
based on the flow quantities between the shock and the wedge, and a length scale L which 
is the distance of interest from the tip of the wedge. Specifically, we introduce coordinates 

A A A A 

Lx, velocities U s u , , time Lt/TJ „ pressure p s U;p> density p s p , temperature T S T and enthalpy 

A. 

h ,h. The governing equations of the flow then become 


where 


dp . x 

— + V • (pu) = 0, 

(2.5a) 

»§= = -Vp + ^[2V • (MS) + £((*' - f /.)■ v • a)], 

(2.56) 

p Dt n w Dt + PrRe y Re 

(2.5c) 

lM]v = pT, h = T, 

(2.5d, e) 

1 f dui duj\ 
e,J 2 dxi ) ’ 

(2.6a) 

$ = 2/j.g : g + in' - jp)(V-u) 2 , 

(2.66) 
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H and fJ,' are the shear and bulk viscosities respectively, both of which have been nondi- 
mensionalized by a typical viscosity jtt g ,P r is the constant Prandtl number, and 

p.U.L 

Re = (2.6c 

P, 

is the Reynolds number which we will assume is large. Henceforth, the subscript s will be 
omitted from M,. We also introduce coordinates (£, t 7 ,js) where £ is measured along the 
shock and 77 is normal to it. Then 

£ = x cos <t> + y sin 0, t] = — xsin<f> + ycos<f>. (2.7) 

For simplicity we will assume (i) that the Prandtl number is unity ( Pr = 1), (ii) that 
the shear viscosity is given by a Chapman law /z = CT, and (iii) that the wedge walls 
are insulating, i.e. = 0 on y — o). At large Reynolds numbers, there is then a simple 
solution to the boundary layer equations in terms of the Dorodnitsyn-Howarth variable 
(e.g. Stewartson 1964), in which the temperature at the wall, T w , is given by 

T w = l + ^(7 — 1)M 2 . (2.8) 

The linear stability of this boundary layer flow has been extensively studied using 
the Orr-Sommerfeld quasi-parallel approximation (e.g. Mack 1969, 1984, 1987). Re- 
cently, Smith (1988) has shown how an asymptotic triple-deck description of lower branch 
Tollmien-Schlichting waves (i.e. ‘first mode’ waves) can be obtained for wave directions 
sufficiently oblique to be outside the local wave-Mach-cone direction. In the limit of large 
Mach number M, Smith (1988) found that the most rapidly growing waves have frequencies 
of order Re*T~ x M ~^ , and wavelengths in the x and z directions of order 

Re~*T£M* and Re~*T£M~* (2.9a, b) 

respectively. As is conventional there are lower, middle, and upper decks in the y direction 
with scales 

Re'sTjM*, Re~sT w and Re~*TjM~* (2.9 c,d,e) 

respectively. 

One aim of this paper is to see how the growth rates of these oblique Tolllmien- 
Schlichting waves are modified by the presence of a shock. To this end, we attempt to 
scale the problem so that the shock occurs in the upper deck. Thus, since <f> must be small, 
it follows from (2. 2d) and (2.9e) that 

<j> ~ etancr ~ Re - * (2.9 /) 
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which from (2.9a) implies an x length scale of order Me tan o. If the complications arising 
from non-parallel effects are to be excluded in order to concentrate on the effects of the 
shock interaction, we require 

Mea « 1, (2.10a) 

i.e. that the wavelength of the Tollmien-Schlichting waves is much less than the distance 
from the apex of the wedge. From (2.2a) and (2.4a), e = 0 (M~ 2 a~ 2 ), so (2.10a) becomes 

oM » 1. (2.106) 

However, for > 0 in (2.4b) we require 

(7 — 1 )M 2 o 2 < 1. (2.10c) 

Hence, in order to consider the effect of the shock on the instability waves without the 
effects of nonparallelism we are forced to make the ‘Newtonian’ assumption 

(7 - 1) « 1. (2.10d) 

It follows that when (7 — 1) is not small there is no systematic asymptotic approach to 

this problem which allows the inclusion of shock effects without those associated with 
non-parallelism. Assuming that (7 — 1 )M 2 is not small, i.e. th^t there is a significant 
temperature variation in the boundary layer, the interaction condition (2.9) becomes using 
(2-8) 

(7 — l^M^a ~ Re*. (2-11) 

Before proceeding to asymptotic expansions based on (2.11), it is convenient to consider 

whether there are any further restrictions on the scales of (7 — 1), M and a. First, we note 
that a pressure/acoustic wave incident on a shock will produce entropy and shear waves 
as well as generating a reflected wave. The entropy /shear waves have a typical lengthscale 
normal to the shock of magnitude 

eo?x~ x ~ Re**^ — 1)*M«<j -1 , (2.12) 

where a is a typical wavenumber in the x-direction (see Section 3 below). For simplicity 
we choose to ignore viscous effects in these waves at leading order, in which case we require 

eaa~ x » (aRe)~?, i.e. from (2.9a) and (2.11) a s M 5 « Re. (2.13) 

If it is assumed that the shock has a viscous internal structure, then its thickness is order 
(ecrRe )- 1 (e.g. Moore 1964), and (2.13) then ensures that the entropy/shear waves have a 
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wavelength much larger than the thickness of the shock. It follows from (2.10a) that the 
waves have typical y scales much less than the width of the upper deck. 

Since nonlinear effects are an important part in transition to turbulence, we will select 
a scaling which leads to a nonlinear problem, before linearizing to obtain an analytic 
solution. Following the scalings in Smith (1988), we conclude that a nonlinear lower-deck 
problem is recovered if 

p ~ (2-14) 

In the upper deck this generates a velocity perturbation normal to the wedge (and normal 
to the shock) of order In order that linearized shock conditions are applicable, 

the undisturbed velocity normal to the shock should be much larger than this, i.e. from 

(2.2e) 

eo»Re~*M~%, i.e. a*M 6 « Re. (2-15) 

If this condition is violated, nonlinear entropy waves result, a difficulty which is not tackled 
here. Also, note that if (2.15) is satisfied, then from (2.10b) so is (2.13). 

In order to fix a scaling we will assume that (cf. (2.10c)) 

(7 — 1 )M 2 a 2 ~ 1, (2.16a) 

then (2.11) and the restrictions (2.10b) and (2.14) imply 

M~o^ReT*, Re~T° « a « Re~**. (2.166, c) 

If the lower-deck is forced to remain linear then the upper bound on a relaxes to o « 
Re~T55. Note also that this scaling implies that the temperature in the undisturbed bound- 
ary is large. For convenience we introduce the scalings 

. . 30 3 8 3 . . 

(7 — 1) = To i Re *, M = oiRei*m, (2.17a, 6) 


then 

T w « -Tmo~ 2 >> 1. 

2 

Also, from (2. 2d), (2.4b) and (2.7) the position of the shock is given by 


y — xe tan a 


^ 3 _ 23 

Re t a t x 


m* 


since 


^ u 2 M 2 


(2.17c) 


(2.18) 
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3. THE TRIPLE-DECK EQUATIONS 


The scalings leading to these equations for compressible flow have been given elsewhere 
(e.g. Stewartson 1974), hence only a brief outline is given here. In all three decks the x, z 
and t scalings are 


x = 1 + aTi2e-r*c!rimU-tX, 3 = a-$Re-?CsTh-*m-*Z, 

4 ^ 1 ^ 1 ^ 1 . 3 

t = a'rRe-'rC'iT w m*\-ST, 


(3.1) 


where A = 0.470 is the Blasius boundary layer skin-friction from the undisturbed middle 
deck solution. 


Lower Deck 

In this layer, the scalings are 

y = a^Re~^C^Twm^Y, u ~ Re~T*C^Twm*\^U, 

v ~ a-TRe-TCiTi\*m~*V, w ~ c-fRe-fCtT^X^W, ( 3 - 2 ) 

p ~ cr~ T' Re~ \ m~ 2 + cr -1 ? Re~$C*\s m~ f P, T ~ T w , p~T~ l . 

On substituting into (2.5) we obtain at leading order 

U X + V Y + W Z = 0, P Y = 0 

U T + UU X + VU Y + WU Z = U YY , (3.3) 

W r + UW X + VW Y + WW Z = -P z + W YY . 

The boundary conditions are 

U = V = W = 0 on 7 = 0, 

(3.4) 

U ->Y + A(X,Z,t), iy — 0 asy-»oo, 

where the conditions as Y — * oo come from matching with the middle deck, and A is the 
so-called displacement function. 
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Middle Deck 


The middle deck has the same thickness as the undisturbed boundary layer, which has 
a finite extent in y when T w » 1 (although an infinite extent in terms of the Dorodnitsyn- 
Howarth variable) . Strictly this means that the middle deck should be divided into three 
regions (i) a boundary-layer region where T » 1, (ii) a region where T ~ 1, and (iii) a 
logarithmically small transition region between the two. 

In region (i) the standard scalings and solutions apply: 

y — Re~*C5T w y*, 


u ~ u* 0 (y*) + Re i*C*T£m * A *Au* 0y ., 
v ~ —<7~ $ Re~ t* C * \ s m~ 5 AxUq, 

(3.5) 

w ~ a-^Re-hc l <\^m-^{R* 0 ul)- l D{X,Z,r), 
p ~ a~ i rRe~^m~ 2 + o~^Re~^C^\^m~^P, 

P~R * 0 {y*) + Re-kAR* 0v ., 

where u^(y) and ii?o(y) are undisturbed velocity and density profiles respectively (note 
that R* 0 = 0(TJ 1 )), and D x = -P z . 

Essentially the same solution holds in regions (ii) and (iii), although minor rescalings 
are necessary. For example in region (ii) Rq = u* 0 = 1, which leads to simplifications in 
the expressions for u and p in particular. Also in region (iii), u* Qy . scales with (logT^)? 
as a result of the kink in the velocity profile when T w » 1. This means that the largest 
velocity perturbations occur in this region. We note that as o approaches R~To this may 
mean that the middle deck becomes nonlinear for smaller wave amplitudes than the lower 
deck. 


Upper Deck 


The scalings for the pressure/acoustic waves here are 


2 ^ 3 3 2. 5 1 

y = a t Re ^C»TwX *m *y, 

(3.6a) 

u ~ 1 -f a t Re 7 C*X 2 m 2 uj, 

(3.66) 

(t/,u/) ~ o~$Re~&C*\5m~5(vi,Wi), 

(3.6c) 

9 



_16_ _ 3 _ 12 _ _4 i,i _ 3 _ 

p ~ a 7 Re 7 m + a 7 Re *p x , 

p ~ 1 + o^Re~^C* 

T ~ 1 + (j$Re~TC*\%m*T\. 

These yield the governing equation for the pressure 

Pixx ~ Piw ~ Pizz 0 . 

One boundary condition comes from matching with the middle deck, i.e. 

Ply = A X x ony = 0 


(3.6d) 

(3.6e) 

(3.6/) 

(3.7) 

(3.8) 


(note p x = P on y = 0), while another boundary condition is to be applied at the shock at 


=Mf) 


_ 3 .5 

2\5 A? 
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C *m 4 


(3.9) 


As is conventional, solutions to (3.7) will be referred to as acoustic waves. 


Shock Conditions 

Wave transmission and reflection across shocks has been studied by Moore (1954), Rib- 
ner (1954), McKenzie and Westphal (1968) and others. In general, whenever an acoustic 
wave is incident on a shock, entropy and vorticity waves are generated in addition to a 
reflected/transmitted acoustic wave. The entropy/ vorticity waves have the same frequency 
and the same wavelengths parallel to the shock as the acoustic waves. However, they prop- 
agate in the direction of the mean flow, which means that for the present scalings that their 
wavelengths normal to the shock are very much less than the acoustic wavelength given by 
(3.6a). The scaling (3.6) does not therefore describe them. In the appendix general jump 
conditions at a shock are given for incident linearized inviscid waves. The limiting process 
appropriate to the above scalings then yields the boundary condition 

Pi = 0 at y = y a . (3.10) 

In deriving (3.10) from the exact solution to the general linear inviscid problem, we ensure 
that proper account is taken of the short wavelength entropy/ vorticity waves. A multiple 
scales approach would be an alternative. 

The entropy /vorticity waves propagate parallel to the direction of the basic flow, which 
is of course parallel to the wedge. It follows that in any situation where the forcing is over a 
distance comparable with the triple-deck scaling (3.1a), i.e. over a distance which although 
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possibly large does not extend asymptotically far upstream, the entropy /vorticity waves 
will be concentrated in a narrow region close to the shock. Hence they cannot affect the 
solution in the middle and lower decks. 2 


4. THE DISPERSION RELATION 

Solutions to the nonlinear system (3.3), (3.4), (3.7-10) can in general only be obtained 
numerically. However, analytic solutions can be found if the waves are of small amplitude 
so that the system linearizes. It is then convenient to focus attention on a single mode, so 
we write 

U~Y + hUe^ aX+ ^ z ~^ + c.c., (V, W, P, A) ~ h{V,W,P, A) e <(«*+**-n*> + c . c> . 

Substituting into (3.3), and linearizing under the assumption that h « 1, we obtain (e.g. 
see Smith, Sykes & Brighton 1977) 


V + PW = 1 4— T Ai{q)dq, 

= iP 2 P f Ai'(f)cfc, 


(4.1c) 


(4.16) 


where 


f = (t'a)3(y — fl/a), f 0 =-* 3 n/a^. 
The solution to (3.7) subject to (3.8) and (3.10) is 


a 2 Binh(( < 9 2 -ot 2 )^(y.-y); 
(/9 2 -a 2 )J cosh((/0 2 -a s )2’i7',) 


P 2 > a 2 , 


a 2 B in((a 2 -/? 2 )*(y.-j7)) ^ g2 < 


(a 2 -^ 2 )* cos((a 2 -/9 2 )2y,) 


Hence, from (4.1), (4.3) and the fact that P = p x on y = 0, it follows that 

’ p 2 > a 2 

(ta) ip a f^At(c)d( _ I ‘ a nh((^ 2 -a 2 )Iv.) 

A ’' ((0) K-ff)*. - fl 2 < a 3 


(a 2 -/9 2 )^ 
tan((a 3 -(9 2 ) ^V.) 


Note that for P 2 > a 2 , the dispersion relation for a flow without a shock is recovered in 
the limit y t — ► oo (Smith 1988). 

2 We note that this is doubly so if a » Re~& since second order viscous effects ensure that the en- 
tropy/vorticity waves decay over a distance much less than the width of the upper deck - although it follows 
from (2.16c) that the lower deck must be linear for our analysis to be valid for such values of a. 
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The dispersion relation (4.4) admits both growing and decaying modes. Neutral curves 
are given by 

Co = (/? 2 - a 2 )* = c 2 «M 2 tanh((/? 2 - a 2 )*y,), (4.5) 


where c 2 « 2.3 and c 2 « 1.0 (a similar expression exists for P 2 < a 2 ). These are plotted 
as solid curves in Figure 2 for three different values of y„. The diagonal dashed curve 
defines the wave-Mach-cone. Above this line the acoustic waves in the upper deck are 
purely sinusoidal, beneath it they either grow or decay in y. When there is no shock in 
the upper deck, waves described by triple-deck theory are constrained to lie in the region 
of parameter space below the diagonal (Smith 1988); with shocks present there is no such 
limitation. 

Asymptotic formulae for these neutral curves can be derived that agree well with the 
numerical results: 


v, \2( n + ?)» y.\ y. ) j 

. (nir \ 2 1 /n7r\ 4 1 1 /n7r\ 2 1 

a ~ /? + — ~ — ) r-53 H — I — ) —w, a> P » 1, n = 1, 2, . . . , (4.66) 

\y s ) 2/3 v y> J S P c *y> \ y» J P~ 


« 1, n = 1,2,. . . , 
(4.6a) 


a 


'J_\ 3 

S*P) ’ 


a << 1, P » 1. 


(4.6c) 


These formulae confirm that there are an infinite number of neutral waves, and that with 
the exception of the subsonic mode, they all asymptote to the line f 3 = a for a, /? large. 

The short dashed lines in Figure 2 also represent waves with zero growth rate, but they 
have an infinite frequency and hence our asymptotic analysis breaks down in their vicinity. 
These lines correspond to the points in the (a, P) plane where 

tan ((a 2 - p 2 )^y s ) -> oo. 


More precisely if we write 


(a 2 - p 2 )*y t = (n+ ^)tt -6 


for n = 0, 1, . . . where 6 is small and positive, then it follows 


n ' = Ij^V^ 2 + ( n + £) 2 fr + • * • » 

(4.7a, 6) 

a- = sgn(<5) (, 3 2 + (n + |) 2 £) ’ \2n R \~* + . . . . 

Hence, if P is held fixed, then when a is within distance 0(6) of a dotted curve, the 
disturbance has frequency 0(6 _1 ) and a growth rate 0(6^). The wave grows or decays 
depending on whether it is below or above the dotted curve respectively. 
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It also follows from (4.7b) that the growth rates in the unstable regions above a — (3 
increase with n. In fact if we hold (3 fixed and let a —* oo, analysis of (4.4) shows 
that the unstable intervals are of range O(oT^) and that the growth rates in the unstable 
intervals are typically 0 ( 0 : 3 ). The frequency in the unstable intervals also increases like as, 
and so the unstable intervals for n » 1 correspond to high frequency, short wavelength 
modes with increasingly high growth rates. However, the fact that these modes occur over 
decreasingly short ranges in a as n inceases means that they might not be excited naturally 
in a physically realistic flow situation. 

In Figure 3 we show the dependence of fl r and fi,- on a for the case y s = 4, (3 = 1. 
We observe the predicted monotonic increase in the growth rates and the shortening of 
the range of unstable wavenumbers as a increases. Results at other values of y t and (3 are 
similar with the range over which the oscillations in a occur decreasing as either y t or (3 
are increased. 

It is perhaps significant that the modes which reduce to those discussed by Smith (1988) 
when y s — > 00 , i.e. the subsonic modes, have growth rates less than those with (3 < a. As 
a measure of the growth rate of the subsonic mode we can take its value close to the line 
a = (3. In fact in all the cases calculated we found that this mode had its maximum value 
in this neighbourhood. In Figure 4, we show the dependence of Q, on y s for different values 
of the spanwise wavenumber (3 with a = .999 (3. We see that increases to a maximum 
and then decreases. Thus, for a given (a, (3) close to the line (3 = a there is an optimum 
value of J7, which maximizes the growth rate. In fact it follows from (4.4) that if 

a 2 — (3 2 = A with |Aj << 1, |/ 4 |A|* << 1, 


then close to a — /? 


.1 T 

y 4 t»a» c* 




Asymptotic analysis for this expression for small and large y t demonstrates that the growth 
rate has a maximum for intermediate y,, e.g. 


fl ~ a 3 y, + (-4r)» as y 4 -> 00 . 

ay, 

We note however that if y,|A|J is not small then the above simplification does not hold. 
Further for y,|A|? large the asymptotic behaviour depends on the sign of A. If A < 0, 
then the eigenrelation (4.4) reduces to 

r = (/? 2 - a 2 ) W(<ro), ( 4 . 8 ) 

J{0 
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which is valid for y t >> 1, a, (3 = 0(1). As expected this is the eigenrelation of Smith 
(1988). In Figure 5 we show the dependence of fl r and fl, on y 3 for different values of 
/ 3 with a — 1. For these subsonic modes the growth rates approach constant values as 
y 3 — ► oo, in agreement with Smith’s (1988) shock-free analysis. Hence subsonic modes of 
fixed wavelength do not feel presence of the shock when y 3 is large (nor do short- wavelength 
subsonic modes when y s = O(l)). 

However, if y t — * oo with A > 0 there is no simplification similar to (4.8) since the 
growth rate oscillates, with the regions of stability and instability becoming increasingly 
thin as y 3 is increased. 

5. CONCLUSIONS 

We have shown that if shock effects are to be included in a Tollmien-Schlichting stability 
analysis of flow past a wedge, then the Newtonian approximation, (7 — 1) << 1, must be 
made if the complicating effects of non-parallelism are to be avoided. With this assumption 
we have seen that the viscous modes have their dispersion relationship crucially altered by 
the presence of a shock at a distance y 3 from the wedge scaled on the upper deck thickness. 
In the limit y 3 — ► 00 the modes of Smith (1988) are recovered and the shock has no zeroth 
order effect on the growth rate of these modes. In addition at large values of y 3 there is an 
infinite discrete spectrum of disturbances which persist to the shock. These modes have 
large growth rates at high frequencies but occur only over small ranges of wavenumber. 
This suggests that the frequency might have to be tuned to produce instability. It is 
therefore an open question as to whether these additional modes play a critical role in the 
transition process in hypersonic flows. We also point out that where our analysis predicts 
infinite frequencies, i.e. close to the dotted curves in Figure 2a,b,c, then our asymptotic 
expansions will fail and a new structure must be set up in order to account for the faster 
disturbance response. This problem has not been investigated in this paper. 

The structure of the eigenfunctions associated with the different modes shown in Fig- 
ure 2 depends on which side of the line a = (3 the given mode occurs. If /3 < a then 
the disturbances are described by a combination of exponentially growing and decaying 
functions. In any of the limits where the shockless eigenrelation is recovered, e.g. y s — ► 00 
or /? — *• 00, the exponentially decaying function dominates and the eigenfunctions also 
tend to those of the shockless problem. The modes which have (3 > a are described by 
trigonometric functions and therefore in any limit involving y s ,(3i or a they are oscillatory 
and 0(1) between the wedge and the shock. This class of the mode is perhaps best thought 
of as sound waves trapped between the wedge and shock and amplified by the boundary 
layer. 
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So far we have only discussed the effect of the shock on viscous modes of instability yet it 
is well-known that compressible shear flows can also support unstable inviscid disturbances. 
The inviscid modes of instability for a compressible boundary layer have wavelengths 
comparable with the boundary layer thickness so we consider perturbations to the flow 
described in §2 but with wavelengths scaled on the main deck thickness. In order to be 
consistent with previous investigations we drop the re-scaling (3.1) introduced in order to 
simplify the triple-deck analysis of §3. 

Following for example Mack (1987), we scale wavelengths on the main deck thickness 
and the wavespeed on the free-stream speed. The pressure perturbation P for an inviscid 
mode then satisfies the compressible Rayleigh equation 


_.. 2u* . T , , 2 a 2 (u — c) 2 M 2 

P" p + =P - (« + p )(l ~ -T^z 7ZZ=-)P = 0. 

u-c T K H A (a 2 + (3 2 )T J 


(5.1) 


Here a dash denotes a derivative with respect to boundary layer variable y*, whilst the 
basic velocity field u and temperature T for an insulated wall are given by 

7-1 


tr = /'(*), T = 1 + M’(l - /»), 


(5.2a, 6) 


where / is the Blasius function and r\ is the Dorodnitsy-Howarth variable defined by 

r « * dy 


H dy * 

,= y. t - 


(5.2 c) 


The quantities a and (3 are downstream and spanwise wavenumbers, whilst c is the 
wavespeed. We confine our attention to modes which satisfy 


P' = 0 on y* = 0, P — ► 0 as y* — > oo, 


(5.3a, ft) 


which together with (5.1) specify an eigenrelation c = c(a, /?). The point at issue here is 
whether the decay of P when y — ► oo is sufficiently rapid to mean that the shock located 
outside the boundary layer has a negligible effect on the disturbance. In order to see 
whether this is the case we discuss the structure of (5.1) in the limit M — > oo. 

Since this work was completed we learned of an independent investigation of (5.1) for 
M » 1 by S. N. Brown and F. T. Smith. They have concentrated on the so-called vorticity 
mode (Mack, 1987), while we will study the acoustic modes. We restrict our discussion 
of (5.1) to the minimum which explains the large M structure of P and shows that in 
the wedge problem considered in this paper the shock generally has only an exponentially 
small effect on the eigenrelation. Further, we concentrate on the generalized inflection 
point neutral modes of (5.1) which have (3 = 0. The discussion we give can be extended 
to unstable two and three-dimensional modes, although an extra critical-layer region then 
needs to be included. 
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The generalized inflection point occurs where u v * tf .T = u y *T y *, and if we write 

where B comes from the large rj asymptotic form of / 

de-in-B )*/ 2 

f — V ~ B -\ — — r= h • . . , 


0 ? - sy 


(6-4) 


(5.5a) 


then we can show that it occurs at 


Y = - 

(7-1 K 

The wavespeed is then given by u evaluated at the inflection point so that (5.1) is not 
singular there. Thus, we have 

c = 1 — ^2 + • • • > (5.56) 

with c = -T-. 

nr— i 

Next we anticipate the change of scale of the boundary layer for M » 1 and write 

A 

CC .n ” 1 " ... » 

M 1 

It follows that the zeroth order approximation to (5.1) for 0(1) values of rj is 
d 2 P 2 f" dP A 2 ^ - 1) 2 (1 - f 2 ) 2 ( 2(1 - f) 


dr) 2 f — 1 dr) 


( l 2ML), = o 

v h-i)(i +n) 


(5.6) 


(5.7) 


For large values of r\ it is found that the solutions of (5.6) are such that 

T) — B + ’ 

where D and E are arbitrary constants, one of which can be fixed as a normalisation 
condition. Anticipating the result we choose D = 1 without loss of generality. E is then 
determined by investigating the solution of (5.1) in the region where Y = 0(l), and also 
in the region above this logarithmically thin layer. In passing it is of interest to note that 
this layer also controls the Gortler vortex mechanism in hypersonic boundary layers (Hall, 
1988). 

For Y — 0(1) we seek a solution of (5.1) by writing 

P = P{Y) + ..., 

and after some manipulation we find that P satisfies 


~ 2 d 2 P ~dP 2 dY 2 dP 

Y 2 ^ + Ys- + = 0, 

dY 2 dY c — dY dY 


(5.8a) 
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(5.86) 


which is to be solved in the range 0 < Y < oo. The solution of this equation is 

~ ~ d 2 Y 2 ~ ~ 

P = D{—— - 2 cdY + c 2 log Y) + E, 
z 

where D and E are constants. This solution is to be matched with (5.7) when Y —> oo and 
a solution of (5.1) valid for r) — B — y/2 log M 2 >> (\/2 log M 2 ) 1 . The former condition 
yields 

n _ 2^2 log M 2 t 

d 2 M 4 + " 

Above this logarithmically thin region there is an outer region with the scaled coordinate 
ff = M~ 2 rj (note from (5.2c) that y* = 0(1) in this region). The solution in this region is 
P ~ E exp(— Arf) = E exp(-ari). Matching this solution to the logarithmic term in (5.8b) 
for Y small we find that 

_ 4c 2 log M 2 

£j . 

ad 2 M 2 ’ 

and matching to the constant yields 

E ~~E. 

Note that this means that the constant term dominates the ^-dependent part of the 
complementary function in (5.8b). However, on substituting (5.8b) into (5.1) we confirm 
that (5.8a) is still the leading order equation for P'(F). An alternative way to see this is 
to work with the equation for the velocity fluctuation normal to the wall, rather than that 
for the pressure fluctuation. 

Finally, matching back to the lower layer adjacent to the wall we conclude that 


E ~ E ~ E « 1. 


Therefore, A is determined by the eigenvalue problem specified by (5.6) together with 

P' = 0 on 77 = 0, P — ► 0 as 77 — ► 00. (5.9) 

This problem was solved numerically to yield the sequence of eigenvalues A = 4.81, 13.84, 
23.34, 33.03, 42.82, 52.67, 62.55, 72.45, 82.37, 92.29, ... . 

In Figure 6 we have shown a comparison of our one term asymptotic form for a with 
the values given by Mack (1987). The agreement is satisfactory except our analysis needs 
modification to explain the ‘kinks’ in the eigencurves found by Mack (but see also the work 
of S.N. Brown and F.T. Smith, private communication). The eigenfunctions associated 
with the above set of eigenvalues are shown in Figure 7. Finally, we note that the higher 
order modes of the eigenvalue problem can be derived by applying the WKB method to 
(5.6). In fact we anticipate that the WKB desription of the modes discussed above coupled 
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to the vorticity mode description of Brown and Smith could be used to explain the kinks 
in the neutral curves given by Mack. Thus we expect that these eigenvalues are split apart 
by an exponentially small amount in the manner discussed by for example DiPrima and 
Hall (1984) for the Taylor problem. 

The main result of the above analysis is that the high Mach number structure of the 
inviscid modes does not lead to a reduced rate of decay at infinity. Hence the shock cannot 
have anything other than an exponentially small effect on them, at least until it is within 
an order one distance of the boundary-layer. Then the steady flow changes and the effect 
of the shock is felt within the outer region. We conclude that the viscous modes are those 
which are most likely to be influenced by the presence of shocks in high Mach number 
flows. 
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APPENDIX: SHOCK CONDITIONS 


In this appendix we give the conditions that must hold at the shock - see Moore (1954), 
Ribner (1954) and McKenzie and Westphal (1968) for similar derivations. 

For a shock given by rj = /(£, z,t), the jump conditions across a shock are (e.g. Majda 
1983) 

I tip] + If (H - M + If IH = °> 

If [H + ff [p « 2 + p] - [pH + If [p*H = o, 

If M + If [pH - [p *> 2 + p] + ff [pH = o, (a.i) 

If [pH + ff [pH - [pvw] + ff [pu > 2 + p] = o, 

If [p£] + ff [«(p£ + p)] - Mp£ + p)] + ffMp£ + p)] = o, 

where £ = I 'i-Tj p + §H 2 + y2 + u; 2 ). The basic solution specific by (2.2) can be shown to 
satisfy these jump conditions with / = 0. 

We assume that there is a small disturbance beneath the shock and write 


(p,u,v,w,p,S) = (R,U,V,W,P,E) + (f,u,v,w,p,e). (A.2a) 

For our flow any waves above the shock propagate towards the shock. Hence the distur- 
bance cannot extend above the shock, where we write 

(p,u,v,w,p,£) = (R,U,V,W ) P,'F). (A.2b) 

The linearized shock conditions are obtained by substituting (A. 2) into (A.l) and neglecting 
all nonlinear disturbance terms. The position of the shock will vary by only a small amount 
from r; = 0, so we write f = f. Also, since the undisturbed flow is a constant above and 
below the shock, to a consistent approximation the jump conditions can be evaluated at 
rj = 0 instead of r? = /. Finally, we assume that the linear disturbance can be expressed 
as a superposition of Fourier modes, so that it is sufficient to study a single mode, i.e. we 
assume 

p <x exp(i(a£ + (3z — fit)), etc.. (A. 3) 

Note that here a, (3 and fi are not the scaled quantites of section 4. After a little manip- 
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ulation it follows from (A.l) that 

-t(fi - aU)(R - R)f + Rv + Vf = 0, 

ia(P - P)f + RVu = 0, 

2RVv + V 2 r + p = 0, 

RVw + i(3(F — P)f = 0, 

-*(ft - aU)(RE - RE)} + iaU(P - P)f + (RE + P)v + VRe + VEr + Vp = 0 
/ and f can be eliminated from the above to obtain 
[3u — aw = 0, 

RVv + p- !2^2Slu = 0, 


(A.4) 


(A.5) 


[^p + -£/>)] s + - £] p + = o, 

where for our basic state and nondimensionalization 


R = l, P = 


7M 2 * 


U = 


(l + e 2 T 2 )2 


r» V = eV, V = 


(1 + e 2 7 2 )*’ 


p _ g g(7-l)(l~e 2 )T 2 - 

7M 2 27(1 + e 2 T 2 ) ’ 

Beneath the shock the linear waves have solutions proportional to e’-i - * 0 *, where k = 
(a,i/,j3) and £ = (£,rj,z). The solutions for the different types of wave have the following 
forms: 


Acoustic: 


(p,r) = (1, *' n< , ( u,v,w ) = 


Entropy: 

Vorticity: 


R(0-U.k) 

n-U.k = 0, r = p=u = v = w — 0 

H -U_.k = 0, (u, u, u;) = (u,, + /9th,), w,)e'-i~ int 


2 _ JP 

R 

(A.6) 

ik.i-int 

w = 0 

(A.7) 

w s )e*i~ int 

(A.8) 


p — r = 0. 

We denote the pressure amplitudes of the incident and reflected acoustic waves by p x and 
p 2 respectively, and their corresponding 77-wavenumbers by *7 and i/ 2 respectively. Then 
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substituting (A.6-8) into (A.5), and using the nondimensional forms of (2.2), we obtain 
after some manipulation 


( Ki + \ _ , ( Ki + K 2 v 2 \ _ _ 

Kn-aV- v x v) Pl + Vn - aU - u 2 v) Pi ~ 


0, 


where 

-(f) - aU){{n - atQ 8 (l + £ 2 T 2 )(3 - 7 + ( 2 + !) g ) - (f! + Philya ± 1) - (1 ± 5 ) g )) 
Kl (O — aU) 2 (l + e 2 T 2 ) + (a 2 + P 2 )eT 2 

_ (7 + 1)(1 -e)eT 
2 (1 + e 2 T 2 )* 

On substituting the scalings (2.16) into the above, we find that at leading order \K\\ » 
\K 2 \, &U » |n — ViV\, and hence 


Pi + P 2 — 0, 


i.e. condition (3.10). Expressions for other quantities can be obtained similarly. For 
a » Re~To the u,w and r perturbations are dominated by the entropy/shear wave 
contributions; in particular for our scalings 


where 


„ 2ctm 2 o 1 r Re** 

U ~ s r Pi, 

(a 2 m 2 -/? 2 )* 

2/?m 2 <7^ Re% 
w ~ — r pi , 

(a 2 m 2 — P 2 ) 2 
—2a m 4 aTiJen 

r ~ — ; — pi, 

(a 2 m 2 — P 2 ) 2 


1_5 _ 3_ „ 

a = a? /te“a, 


23 3 ~ 

P — a 7 Rtf p. 
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The geometry of the wedge and shock for a high Mach number flow. 

The neutral curve a = a((3) for y g = 1. 

The neutral curve a = a (/?) for y t — 4. 

The neutral curve a — a(j3 ) for y t = 16. 

The growth rate fi,- as a function of a for f/, = 4,/? = 1. 

The wavespeed fi r as a function of a for y, = 4,/? = 1. 

fl r and 0* as functions of y t for = 1, a = 0.999/?. 

Q r and fi, as functions of y t for (3 = 2, a = 0.999/?. 

The growth rates as a function of y g for a = 1,/? = 2,3, 4. 

The neutral curves for the generalized inflection point modes with /? = 0. 
Mack’s (1987) results. High M asymptotic prediction. 

The eigenfunction P(y) for A = 4.81. 

The eigenfunction P{y) for A = 13.84. 

The eigenfunction P(y) for A = 23.34. 

The eigenfunction P(y) for A = 33.03. 

The eigenfunction P{y) for A — 42.82. 

The eigenfunction P(y) for A = 52.67. 

The eigenfunction P{y) for A = 62.55. 

The eigenfunction P(y) for A = 72.45. 
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